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ABSTRACT 


Energy distribution relations which describe the desorption behavior from 
heterogeneous surfaces are presented. A graphical method of obtaining these is de- 
scribed, and the results are compared with a more exact computer -determined solution. 
Adsorption isotherms representing equilibrium coverage for the systems H 2 ~W, CO-W, 
N 2 ~W, CO-Mo, and N 2 ~Mo are presented covering a pressure range from 10" 12 to 
1 torr, a fractional surface coverage from 0. 02 to 0. 98 and for the appropriate temper- 
ature range in each case. 



I 


ADSORPTION-DESORPTION BEHAVIOR OF HOMOGENEOUS 

AND HETEROGENEOUS METAL SURFACES 
by Thai neW. Reynolds 

Lewis Research Center 
SUMMARY 

Equations and curves describing the adsorption of gas on metal substrates under 
equilibrium conditions, and the desorption of the gas from the surface under a linear 
temperature increase of the surface temperature are presented. 

Results of an experimental investigation using the flash-filament desorption tech- 
nique are used to obtain an energy -distribution relation for several gas -metal surface 
combinations. A graphical method of approximating the energy distribution relation 
which describes the desorption behavior from heterogeneous surfaces is presented and 
compared with the more exact solutions obtained on a digital computer. The range of the 
desorption energy spread was about 30 to 40 percent either side of the desorption energy 
value giving maximum desorption rate. 

Systems for which quantitative experimental data are presented are hydrogen- 
tungsten, carbon monoxide -tungsten, nitrogen-tungsten, hydrogen-molybdenum, carbon 
monoxide -molybdenum, nitrogen-molybdenum, and hydrogen-tantalum. Some observa- 
tions on the behavior of other gas -metal combinations are also made. 

The energy distribution relation is then used to calculate adsorption isotherms. 

These isotherm plots express the equilibria expected over a range of pressures from 
-12 

10 to 1 torr, for fractional surface coverages from 0. 02 to 0. 98 monolayers and for 
the appropriate temperature range in each case. 


INTRODUCTION 

Study of the adsorption and desorption behavior of gases on surfaces has become of 
considerable interest in recent years because of the availability and extended use of 
high vacuum equipment. In fact, the attainment of ultrahigh vacuum is usually accom- 
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plished by taking advantage of surface adsorption phenomena, such as in cryosorption 
and sublimation pumping. 

Two simple M rules-of -thumb", which are frequently cited, serve to illustrate the 
relative importance of surface adsorbed species to the attainment of high vacua; they 
are the following: 

(1) There may be 3x10 times as many particles on 1 square centimeter of an ad- 

O 

sorbed monolayer as there are in 1 cubic centimeter of the gas phase at 10 ~ torr. 

(2) Gas molecules arrive at a surface at a rate equivalent to about one monolayer 
per second at 10 torr. 

Consideration of the first point demonstrates that very small changes in surface 
coverage can involve enough particles to affect the pressure level of a system con- 
siderably. 

The second point illustrates the necessity for working at very low pressures when 
studying "clean surface" phenomena. 

In addition to the effect of adsorbed gas on vacuum system behavior, there are many 
surface phenomena that are greatly affected by the presence of adsorbed gas layers of a 
monolayer or less. Electron emission, surface ionization, surface diffusion, and 
catalysis are a few such phenomena. 

One important experimental technique for studying surface adsorption and desorption 
behavior is that of flash -filament desorption (ref. 1). This technique involves the 
thermal desorption of adsorbed gases by heating the surface to high temperatures. How- 
ever, work with the field emission microscope has pointed out that considerable surface 
migration of metals occurs at temperatures of one -third the melting point and lower 
(ref. 2). Thus, unless a surface which has been heated above this temperature level 
tends to reform to some equilibrium structure upon cooling, one would not expect to be 
investigating one particular surface structure. The temperatures required to clean a 
surface and to flash off chemisorbed material are, in general, considerably higher than 
the temperatures at which the surface may become quite mobile. 

Considerable literature data are available for the flash filament desorption tech- 
nique (ref. 1). The reproducibility of these data after successive flashes to high tem- 
perature levels must constitute some evidence for the return to an equilibrium surface 
structure on cooling. (Equilibrium, in this case, means only that the surface attains a 
structure which has repeatable desorption behavior.) 

This report discusses adsorption isotherm relations and transient desorption be- 
havior of energetically homogeneous and heterogeneous surfaces in general. It discusses 
the determination of site -energy relations for heterogeneous surfaces from experimental 
flash-filament desorption curves, and finally, presents some calculated adsorption 
isotherms which are based on the flash-filament-determined surface desorption energies. 
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The term, isotherm, is used herein to designate the relation between fractional 
surface coverage, surface temperature, and arrival rate or pressure. 


THEORY 

Equilibrium Adsorption- Desorption Isotherm Relations 

Consider a surface immersed in a gas and at equilibrium with respect to adsorption 
and desorption of the gas. 

The desorption rate of a gas from a constant energy surface is usually described by 
the relation (refs. 3 and 4). 

_n -E/kT 

f = — — e 8 (1) 


where 

2 

v desorption rate, particles/(cm )(sec) 

2 

o surface coverage, particles/cm 

n order of desorption relation 

"r „ constant, sec 

on 

E desorption energy barrier, eV 

k Boltzmann constant, eV/k 

T g surface temperature, K 

(Symbols are defined in the appendix. ) 

The adsorption rate is generally expressed by defining a sticking coefficient S such 

that 


Adsorption rate = jiS 


( 2 ) 


where [i is the arrival rate in particles per square centimeter per second. 

At equilibrium the adsorption and desorption rates are equal, so that a simple iso- 
therm relation becomes 


u n e " E / kT s 
T on 


(3) 
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A more generally useful form for isotherm information has the coverage expressed 
explicitly as 


o = tin, T s ) (4) 

Then, values for S, r Qn , n, and E must be known in order to use equation (3) to des- 
cribe the adsorption data in the form of equation (4) . 

While the surface may be energetically uniform, E can vary with coverage because 
of interaction among the adsorbed particles. In the analytical treatment of uniform 
surfaces herein it will be assumed, as is usually done(e. g. , ref. 3), that E does not 
vary with coverage. The value of r 1 is sometimes treated as equal to h/kT or 
sometimes as the constant, 10“ second. Also values of and E may both be deter- 
mined from plots of the experimental desorption data over a range of temperature. Each 
method may yield a different value of E. A discussion of the significance of the term, 
r , may be found in reference 4. 

Experimental values of E may be deduced readily from flash-filament desorption 
data as will be discussed shortly, but it should be kept in mind that the values reported 
depend on the equation used to obtain them. Desorption energy, per se, is not measured 
directly. 

The reaction order n expresses the dependence of the desorption rate on the sur- 
face concentration of adsorbed molecules. The desorption rate of a monatomic gas or of 
a molecular gas which does not dissociate upon adsorption would be expected to depend 
directly on the surface concentration (i.e. , n = 1). For a molecular gas which dissoci- 
ates upon adsorption and recombines during desorption, the desorption rate would depend 
upon a higher order of surface concentration (n > 1). For most adsorbed species of 
interest, the reaction order n will be either 1 or 2 - or possibly some intermediate 
value if conditions are such that more than one desorption reaction is occurring simulta- 
neously. 

The sticking coefficient S is the remaining factor to be evaluated. It may vary with 
both temperature and surface coverage. Experimental values of S are required since 
there is no theoretical basis for calculating a value of sticking coefficient. Further- 
more, there is considerable variation in reported values of sticking coefficient for the 
same systems. 

For situations in which the temperature of the surface is above that corresponding 
to the saturation temperature of the gas, the adsorption amount is frequently found to 
approach a limiting coverage of a monolayer (ref. 5). One way of incorporating this 
limit into an isotherm relation is to assume that the sticking coefficient variation with 
coverage is simply 
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(5) 


S = S Q (1 - 0) 

where 6 = ct/ct^ and cr m is the value for monolayer coverage. This is equivalent to 
assuming that particles arriving at filled adsorption sites are reflected. 

Using this relation results in the isotherm equation (3) becoming 


M( i - s)s 0 


c^ e - E / kT s 


( 6 ) 


This can be rearranged to yield the forms 


e 

l - e 



eAt s 


for n = 1 


(7) 


and 



t o2 E / kT s 

= -2£ e s 

^ o 2 
a„ 
m 


for n = 2 


(8) 


Theoretical isotherm curves based on equations (7) and (8) are shown in figure 1 for 

a range of variables. Figure 1 has been made in the form of a nomograph for more 

13 

compact presentation. In constructing these curves, values of t ^ = 10 second and 
t Q 2 = 50 seconds per square centimeter were used. The value of here is taken to 
be consistent with the ratio t o2 / t o1 = CT m / 2 determined from reference 6 and cr m is 
taken as l.OxlO 15 particles per square centimeter. 

A nomograph relating arrival rate, pressure, temperature, and molecular weight 
is presented as figure 2 so that arrival rates needed for figure 1 are readily obtained 
for any particular gas. The relation used for figure 2 is 

p = — 54 xl C> P. molecules/(cm^)(sec) (9) 

VT g m 

where p is the pressure in torr and m is the molecular weight. The combination, 
p/p, is used later in the report, and it should be noted that this expression is independ- 
ent of pressure. 
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Transient Desorption Relations 


Homogeneous surfac es. - One technique that has proved very useful for investigating 
the kinetics of desorption is the ’’flash -filament” method. The theory associated with 
this technique has been adequately presented in references 3, 7, and 8. The method 
consists of heating a surface, upon which a gas has been adsorbed, at a programmed 
rate, and noting the desorption rate as a function of temperature. Various time- 
temperature schedules are discussed in the references, but only one, a linear rate of 
temperature rise, will be considered here. 

The desorption rate equation for the uniform surface is 


n -E/kT 
v = _ e s 

T 

on 


(i) 


Since the desorbing particles come only from the adsorbed layer 


v - - 


da 

dt 


( 10 ) 


and equation (1) becomes 


^ -E/kT £ 


da_ cr~ — s 
dt r 

on 


(ID 


If a linear rate of temperature rise 


P = 


dT 
£ 

dt 


( 12 ) 


is assumed, equation (11) can be expressed as 


da 1 - E / kT s 
- - 1 e s dT 

o" 0 T on 


( 13 ) 


In order to integrate equation (13) the following approximation is found to be useful: 
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d 


( 14 ) 


dT 

s 




E * E AT b 

-■ V 

k 



2 

E/kT g 


For the present situation, only values of E/kT 20 will be of much interest, and 

s 

the approximation that 2/(E/kT g ) is negligible compared to unity is made. With this 
approximation, equation (13) can be integrated to 


and 



a ^ T ol 


, / 9 -E/kT 
k/ T 2 e s 

El s 




for n = 1 


(15) 


1 

cr 


J_ _ 1 k/ T 2 e " E/kT ; 
a 0 / 3 t o2 E \ S 




for n = 2 


(16) 


As long as T g ^:l.lT g Q, the T g0 term will be negligible with respect to the T g term 
in equations (15) and (16). Dropping the T g Q terms and rearranging yield 


— = exp 
a 0 


! kT l -EAT g 
^ T ol E 


n = 1 


(17) 


and 


_o_ _ 1 

a 0 ct 0 kT g -E/kT 

i + — y e 

^ t o2 E 


n = 2 


(18) 


Equations (17) and (18) are shown plotted in figure 3 for a range of variables E 
and T and for several values of j3. Equations (17) and (18) describe the change in 

o 

surface coverage, expressed as the ratio of coverage at any time (or temperature) to 
initial surface coverage ct/ctq for a surface whose temperature is increased at a con- 
stant rate. There is in general a narrow temperature range over which most of the 
desorption occurs for a constant energy surface. This range does not vary too greatly 
even for large changes in £, so that it is only necessary that the rate of temperature 
change of the surface be relatively constant over the portion of the temperature band 
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through which desorption occurs, to be able to interpret experimental desorption in terms 
of these theoretical curves. 

It may be noted from equation (1) that, as temperature rises, the desorption rate 
increases. As desorption occurs, the adsorbed layer becomes depleted so that a maxi- 
mum desorption rate finally occurs, after which the desorption rate drops as the tem- 
perature continues to rise. The desorption rate for the linear increase in temperature 
shows the typical trend illustrated in figure 4. 

The theoretical temperature T max at which the maximum desorption rate occurs 
is shown in figure 5 for both first- and second-order desorption, for a range of desorp- 
tion energy values. 

Heterogeneous surfaces . - Analysis of adsorption -desorption behavior for energet- 
ically heterogeneous surfaces - that is, surfaces for which E varies over the area 
independent of any variation that may be due to interaction between adsorbed particles - 
has been treated in references 3 and 7. 

The desorption rate from such a surface is obtained by summing up the area- 
weighted rates of the separate patches with individual desorption energy E. If the 
distribution of energetic sites can be considered a continuum, then a site-energy distri- 
bution function f(E) can be used and the coverage ratio <r(T )/a n for a heterogeneous 

s u 

surface can be expressed as 



These relations cannot generally be evaluated except for certain simple situations. 
Numerical methods of deducing the distribution function from experimental data and 
equations (19) and (20) must be employed as discussed in reference 7. A simple graphi- 
cal method for approximating the distribution function is discussed in the next section of 
this report. A comparison of the results of this graphical method and the numerical 
solution is presented. Of course, the same restriction must be applied to the site- 
energy distribution relation as was put upon the desorption energy itself, namely, that 
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its values will depend upon the relations assumed to calculate it. The distribution may 
vary from one sample to another of the same metal. 

Site-Energy Distribution 

When a desorption trace for an initial coverage of a monolayer is available, a site- 
energy distribution characterizing the desorption can be determined. Assume a typical 
plot of experimental or/cr m against temperature superimposed on the set of desorption 
curves for a homogeneous surface as in sketch (a). 



(a) 


If the experimental range of temperatures and energies encompassing the complete 
desorption curve is broken up into a number of segments, a numerical solution to equa- 
tion (19) or (20) can be obtained by writing a set of simultaneous linear equations. 
Equation (19) or (20), for instance, becomes 


E=E, 


V f(E) — 

a—/ cr 




AE = 


ffCrj 


( 21 ) 


E=E, 


m 


m 


The right-hand terms of equation (21) are obtained from the experimental data curve. 
The ct(E, T_)/cr relation used is either equation (17) for n = 1 or equation (18) for 
n = 2. 

A difficulty in the attempted solution of this set of simultaneous equations arose 
because of the exponential nature of the cr(E, TJ/cr terms. Negative values of some 
of the f(E) terms would arise in order to satisfy the solution. These are, of course, 
physically meaningless in this instance. Herein, an iterative method of solution was 
used. The iterative process involved correcting the value of f(E) for the E which 
gave the maximum desorption rate at the particular temperature. In this iterative 
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solution method, no negative f(E) values were permitted to occur. However, to obtain 
a reasonably smooth continuous solution may require use of a large number of segments. 
For the solutions presented herein, a 36 x 36 matrix of points was used. The resulting 
energy -distribution relation reproduced the experimental data curve to well within the 
accuracy of its initial determination. 

It is also of interest to note that an approximation to the energy -distribution curve 
may be obtained directly from the graph of sketch (a). The method of doing so will be 
outlined here, and a comparison of this graphical method with the previous iterative 
solution will be shown in the RESULTS AND DISCUSSION section of this report. 

Since the experimental curve of sketch (a) has been obtained for a monolayer cover- 
age, the area fraction can be made to correspond with the coverage ratio. Furthermore, 
the relatively steep constant energy desorption curves give credence to the concept of the 
adsorption sites filling up and desorbing sequentially. It can be supposed then, that, 
since (in sketch (a)) no adsorption has occurred at about 500 K, there are no adsorption 
energy sites less than 1. 5 electron volts. Similarly, since desorption is essentially 
complete by 1600 K, there are no sites with interaction energy greater than about 
3. 5 electron volts. Assuming the sequential filling up of sites, then, one can determine a 
value of F(E), the fraction of area which has a desorption energy equal or greater than 
any selected value E, by noting the value of E corresponding to <r/a (in this case 
assumed equivalent to F(E)) along the experimental curve. A continuous curve of this 
area function against energy can be drawn (see sketch (b)). 



(b) 


The site-energy distribution function f(E) and the function F(E) described pre- 
viously are related by 


f(E) = - 

dE 


( 22 ) 
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It might be helpful to examine the desorption curves obtained from some assumed 
site -energy distribution functions. Two such distributions are assumed: (I) - two equal 
areas of only two different energy sites, E = 2. 5 and 4. 5 electron volts, and (n) - three 
equal areas of just three energy sites, E = 2. 5, 3. 5, and 4. 5 electron volts. Sketch (c) 
shows the appearance of the F(E) against E plots for these assumed distributions. In 
sketch (d) the corresponding comparative cx/o^ against temperature curves that would 
be obtained from these distributions is shown. 



The trend is evident; as the fractional areas of the patches of any particular energy 
level decrease, the inflections in the <x/(r against temperature curve becomes less 
pronounced. 


Isotherm Determination 


Once a site -energy distribution relation for a particular gas -surface combination 
has been determined, the adsorption isotherms are readily calculable. If the distribu- 
tion function were expressible analytically, the coverage fraction 0 would be obtained 
as 


°<T S ) = 


jf 


f(E)0(E, T ) dE 
0 s 



OO 

f (E) dE 


(23) 


Herein, the total coverage fraction was evaluated numerically using the 36 -point ap- 
proximation to the relation f(E) 
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( 24 ) 


36 

£ 

0(T S ) = hi 


f(E.)0(E.,T s ) AE i 


E f (E i } AE i 

i=l 


The values of 0(E, T g ) are obtained from equations (7) or (8) for n = 1 or 2, respec- 
tively. 

EXPERIMENTAL APPARATUS AND PROCEDURE 

A schematic layout of the apparatus used to obtain the experimental data reported 
herein is shown in sketch (e). 


To 



gage 

(e) 

The test surfaces were wires 0. 01 inch (2. 54x10"^ m) or 0. 005 inch (1.27x10”^ m) 
in diameter by 16 inches (0. 406 m) long mounted inside a test volume of 1047 cubic cen- 
timeters. A T -valve between this chamber and the vacuum system permitted control of 
the conductance or pumping rate from the chamber. A Bayard -Alpert type ionization 
gage, operated at 0. 1 normal emission level to reduce gage pumping effects, was used 
to measure test chamber pressure. Gas flow could be introduced to the chamber at 
controlled rates through a calibrated porous plug. A mass spectrometer was used to 
identify background constituents and desorbing species. All parts of the system were 
equipped with bakeout heaters . 

The test filament was heated by a direct-current electric supply, which in turn was 
controlled by a motor driven transformer in order to achieve the desired temperature 
rise rate. The filament voltage and current readings, used to determine the wire tem- 
perature, were monitored continuously, along with ionization gage and mass spectrom- 
eter outputs, on a recording oscillograph. 
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Calibrations 


The various volumes of the rig were determined. The porous plug was calibrated by 
measuring the rate of decrease of pressure in a known volume upstream of the plug. A 
typical plot is shown in sketch (f ) for carbon dioxide (COg). 



<f) 


The plug was fabricated from 2 to 3 -micron tungsten (W) powder so that the pore 
sizes were of this order of magnitude. The flow calibrations were made at pressures 
below 20 millimeter upstream pressure, at which level flow through this plug would be 
expected to be entirely within the free-molecule range (ref. 9). 

The time constant r was calculated from the relation 


£2 = g-At/r (25) 

Pi 


For free-molecule flow, the constant for a system can be derived as 


T = 


2. 73X10" 4 



sec 


(26) 


and is seen to be proportional to V™"- 

A plot of the time constants obtained for nitrogen (Ng), hydrogen (Hg), and helium 
(He) confirm the -\/m dependence for this plug in the range used (sketch (g)). 
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(g) 


The conductance of the T-valve which is needed to obtain quantitative desorption 
rates, was determined in terms of the ion gage readings for several gases. A flow rate 
was set through the porous plug, and the ion gage was then read for various settings of 
the T-valve. This was repeated for several flow rates. A typical plot of these readings 
for Ng is shown in sketch (h). 


8xlCr 5 



(h) 


From curves such as these then, the flow rate through the T-valve can be deter- 
mined for any setting of the valve and reading of the ion gage. It should be noted that 
the absolute gage factor of the ion gage, which would be needed to obtain true system 
pressures, need not be known to obtain true flow rates in this case. Only if one now 
wishes to determine flow rate through the T-valve for some gas for which a primary 
calibration was not obtained will it be necessary to know the relative gage factor for that 
gas compared to one used in the calibrations. 
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The wire temperature was determined from its resistance using the resistivity data 
of reference 10 for tungsten and reference 11 for molybdenum (Mo) and tantalum (Ta). 
Temperatures in the higher temperature range were checked initially using an optical 
pyrometer and also by comparing experimental power requirements against calculated 
radiant energy loss from the wire. All methods were found to be in excellent agreement. 
A calculation of the expected temperature distribution over the length of the wire was 
made following the relations noted in reference 12 . The cooling effects on the wire ends 
caused by conduction to the supports reduces the area which is uniform in temperature 
(i. e. , to within experimental accuracy) to about 95 percent of the wire length. 


Procedure 

A typical procedure to obtain a desorption trace was as follows. The wire was 
flashed to clean its surface, after which it was returned to room temperature. Gas pres- 
sure was established at some specified pressure level and adsorption allowed to occur 
for a certain length of time, herein called "adsorption time". The T-valve was then 
opened to a setting for which the transmission area had been determined. During this 
time, the system pressure would become established at a value designated p . 

CC[ 

The temperature of the wire was then increased at a specified temperature rise rate 
(3 to a temperature level at which the particular gas would be all desorbed. The pres- 
sure and temperature curves as functions of time were followed during this desorption 
step. Simultaneous mass traces were also taken during desorption to identify the 
desorbing gas. 

These data were then analyzed in the manner outlined in the next section. 


Data Analysis 

Consider the system shown schematically in sketch (i) which incorporates the 
essential features of the experimental apparatus. 



(i) 


15 




A mass flow balance on the system yields 


dN 

dt 


F c- F d + 


F fc - F 
d a 


(27) 


where 

N total number of atoms in chamber 
F c desorption rate from test surface, particles/sec 
F d adsorption rate onto test surface, particles/sec 

F^ combined flow rate into chamber from all sources other than test surface (i.e. , 
leak, outgassing, intentionally introduced flow, etc.), particles/sec 

F_ net flow rate out of chamber to pump, particles/sec 
a 

These various flows are determined as follows: 

F a = p C A<P c " Pexh’ (28) 

where is the equivalent area of a sharp-edged orifice which has the same conduct- 
ance as the system has through the opening to the vacuum pump. The rate F_ will 

a 

either be set at some predetermined value or be known from calibration procedures. 
Then 


F c = vA g (29) 

and 

F d = mSA s (30) 

where A„ is the surface area, herein used as irdL where d is wire diameter and L 
s 

is wire length. Substituting equations (28), (29), and (30) into equation (27) yields 

f^A s -MSA s+ F b -ttC A (p c -p exh ) (31) 

By using the ideal gas law, equation (32) is put in terms of pressure as 

f = [M s - pSA s + F b - t C A (p c - p exh )j (32) 
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This relation (eq. (32)) describes the behavior of pressure in the test volume under 
the combined actions of adsorption and desorption from the test surface, introduced gas 
flows, and continuous pumpout of gas from the chamber. Interpretation of the system 
behavior for several different ways of operating will be apparent using equation (32). 

In determining the transmission area of the T -valve for instance, a steady flow 
rate of gas was introduced. This rate is much greater than any leak or outgassing rate, 
so that is equal to the known introduced rate in this case. The pressure reached an 
equilibrium value p so that dp/dt = 0. The surfaces reach equilibrium with respect 
to adsorption and desorption rates; that is, 

vA s = MSA g (33) 

Thus, 



Now consider the desorption trace in sketch (j), which might be typical of the pres- 
sure-time behavior in the system under continuous temperature increase of the surface, 
with continuous pumpout. 



U) 


Since in equation (32) represents the combined flows of introduced gas and 
background leaks and outgassing rates, equation (34) applies whenever the surface is at 
equilibrium with respect to adsorption and desorption rates and dp/dt = 0. Substitution 
of equation (34) into equation (32) yields 
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(36) 


toe 

dt 


kT 


g 


vk B - " SA s - C A ; (Pc - Peq) 


The rate at which particles are removed from the system over and above the rate 
as a result of continuous temperature increase of the test surface is given by 


dN 

dt 




(37) 


Since it is assumed that these particles only come from the test surface and that read- 
sorption is negligible during the desorption phase, the total integral under the pressure - 
time curve above p e ^ yields the total initial surface coverage 


Vs 




C A (p c - p eq> dt = 


c 

h 


^ C A ( Pc-Peq) dt 


(38) 


At some intermediate time, the number of particles that have been desorbed involves 
the area under the curve up to that point in time, plus a contribution which consists of 
the particles that have been desorbed but not yet pumped out. This latter term is 


AN. 


desorbed but 
not pumped 


= — (Pn " PoJ 


kT 


g 


eq' 


(39) 


Thus at time t, 


a(t)A s = <t 0 A s - fv A s dt 

•An 


(40) 


Using equation (38) results in 

*t„ 


f 

r l 


a(t)A s = I ; C A (p - p eq )dt 

r 


/■ 

J t, p 


- C A (P - P eq> dt - <P - W (41) 


The pressure -time curve during desorption may be analyzed by means of equa- 
tion (41) to yield a surface concentration-temperature relation, which may then be com- 
pared with the theoretical curves of figure 3. If the pumping rate through the T -valve is 
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sufficiently rapid in comparison to the desorption rate, then the temperature at which the 
pressure reaches a maximum corresponds closely to the temperature at which the maxi- 
mum desorption rate occurs (ref. 3). For the conditions of this experiment, the maxi- 
mum difference between these two temperatures is about 2 percent according to the 
analysis of reference 3. 


RESULTS AND DISCUSSION 

Results are presented for each particular metal surface investigated and for each 
gas for which data were obtained with that metal. Site -energy distributions calculated 
from these data are presented. Finally, some calculated isotherms are presented using 
these site -energy relations. 


Adsorption- Desorption Data 

Generally after bakeout of the system and heating of the wire filaments to a high 
temperature level (2800 K for tungsten and 2500 K for molybdenum and tantalum), the 
background gas composition consisted mainly of Hg, He, water (^O), and carbon 
monoxide (CO). In no instance in the tests herein was any He or HgO desorption noted 
from the test surfaces. 

A series of curves, typical of those obtained with most systems, will be shown 
only for the first gas -surface combination presented. Data for the remainder of the 
systems will be summarized in table I or in the figures. 

Tungsten . - Tungsten is the easiest metal on which to obtain a clean surface since 
it can be heated to a higher temperature than any other metal without melting or sublim- 
ing. 

Hydrogen: Desorption of hydrogen from tungsten is shown in sketch (k) for a H 2 
background pressure of 7. 8x10”^ torr during adsorption. The area under these curves 
is proportional to the initial surface coverage. The surface coverage buildup that oc- 
curred is shown in sketch (Z ) . 

The maximum desorption rate is seen (in sketch (k)) to occur at lower temperatures 
as the initial coverage increases. This behavior is typical of a second-order desorption 
for constant E surfaces as is illustrated by the theoretical curves of figure 5(b). A 
plot of the values of this maximum temperature T max for this data and for another 
similar set of Hg-W data is shown in figure 5(b). The general level of the desorption 
energy is seen to be 1. 7 to 2. 0 electron volts assuming a second-order desorption. 

Table II contains a summary of some literature values (ref. 1) with which results on the 
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various systems reported herein may be compared. 

A typical desorption curve obtained after adsorption at hydrogen pressures up to 

10 -5 torr is shown in figure 6 along with other typical curves. The maximum surface 

15 

coverage obtained in these tests was 1.0x10 H-atoms per square centimeter. This 
value is taken as monolayer coverage for the H 2 ~W system. The data curve of figure 6 
was used to determine the site -energy distribution relation for the Hg-W system in the 
manner previously described. 

Carbon monoxide: The maximum amount of adsorbed CO on tungsten measured was 

5. 3x10^ CO-molecules per square centimeter. Figure 6 shows the desorption trace for 

this maximum coverage case. Since CO apparently does not dissociate on adsorption 

(ref. 13), the temperatures of maximum desorption rate for this system are plotted on 

the first-order theoretical plot (fig. 5(a)). The general level of desorption energy for 

the CO-W system is seen to be 3. 3 to 3. 4 electron volts. 

Carbon monoxide desorption curves obtained both from adsorption of background 

CO gas and from introduction of C0 2 into the system which dissociated on adsorption. 

There was no difference in desorption characteristics of the CO between the two cases. 

Nitrogen: Desorption of nitrogen from tungsten is also shown in figure 6 for the 

14 

maximum surface coverage of 6. 44x10 N-atoms per square centimeter. The temper- 
atures of the maximum desorption rate obtained for the N 2 -W system are plotted in 
figure 5(b). The desorption energy level is 3. 7 to 4. 1 electron volts. 

Other gases: Carbon dioxide and oxygen were also introduced into the system but 
no quantitative adsorption data were obtained with these gases on tungsten. Qualitatively, 
the following behavior was observed. It appears that C0 2 is dissociated readily on 
contact with the tungsten surface, even at low temperatures (i.e. , near ambient). This 
has been discussed elsewhere (ref. 13). A small amount of C0 2 desorbed as such at 
temperatures around 500 K. The remaining adsorbed material desorbed as CO. The 
excess oxygen may either have reacted with carbon already in the metal, coming off as 
CO, or remained to desorb as O or oxygen (0 2 ) at higher temperature levels. Similarly, 
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no quantitative data were obtained on C> 2 desorption. The only qualitative reproducible 
behavior noted was a desorption peak which coincided with the 32 mass line and occurred 
in the temperature range 2340 to 2440 K. If the kinetics for desorption are first order, 
a desorption energy of 6 to 6. 3 electron volts is indicated (fig. 5(a)). If the desorption 
is second order, the desorption energy could range down to around 5. 5 electron volts 
(fig. 5(b)). Since coverage was not determined, the actual value of ct q //3 could not be 
determined. 

Molybdenum . - Molybdenum cannot be heated to quite as high a temperature as 
tungsten. Nevertheless, it can still be cleaned of most species by high temperature 
treatment. It has been reported (ref. 14) that CO adsorbed on molybdenum cannot be 
removed by high temperature treatment below the melting point of molybdenum. Data 
presented in the following sections are in contradiction to this statement. 

Hydrogen: Temperatures at which the maximum desorption rate occurred with the 
H 2 -Mo system are shown plotted in figure 5(b). If second-order desorption is assumed, 
the energy level is 1. 35 to 1. 55 electron volts. 

The maximum amount of adsorbed H 2 measured was 2. 65X10 1 " 1 H-atoms per square 
centimeter. This maximum measured amount here is undoubtedly not monolayer cover- 
age, but more probably was close to the equilibrium value expected in the system at 
the existing H 2 pressure and temperatures. Without being able to take the surface tem- 
perature below ambient in this experimental setup, it is impossible to reach and 
maintain monolayer coverage. Even though higher coverages may be reached by raising 
the H 2 pressure, the areas of lower desorption energy will rapidly desorb during the 
time interval in which the system pressure is being lowered prior to starting a desorp- 
tion trace. So, while the adsorption of hydrogen on molybdenum could be followed over 
some of the area, desorption from a complete monolayer could not be observed in this 
experiment. 

Carbon monoxide: The temperature at which the maximum desorption rate occurred 
for the CO-Mo system did not vary with the initial adsorbed amount. This is typical of 
the expected behavior of first-order desorption from monoenergetic surfaces. The 
values of T are consequently shown on the theoretical first-order graph (fig. 5(a)) 

and show the desorption energy level to be about 3.1 to 3.3 electron volts. 

14 

The desorption trace for the maximum coverage case, 4. 83x10 CO-molecules per 
square centimeter, is shown in figure 6. This number is treated as monolayer concen- 
tration herein for the CO-Mo system. 

Nitrogen: A desorption trace for the maximum coverage case for nitrogen on 

15 

molybdenum is shown in figure 6. The coverage, 1. 16x10 N-atoms per square 
centimeter, is taken as the monolayer value for this system. A plot of the temperature 
for maximum desorption rates, shown in figure 5(b), yields the general level of desorp- 
tion energy of 3.0 electron volts for N 2 -Mo. 
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Mercury: Several attempts to measure adsorption of mercury on molybdenum 
yielded no indication of adsorption in the range above room temperature. This was sur- 
prising in view of the results reported in reference 15, which indicate adsorbed amounts 
at temperatures up to 650 K, a readily measurable condition in the present experiment. 
The fact that no measurable adsorption was noted at room temperature would indicate a 
range of desorption energies for the Hg-Mo system below about 1. 0 electron volt (see 
fig. 3). Reference 15 cites a desorption energy value of 1. 8 electron volts. 

Tantalum . - Tantalum exhibited adsorption behavior quite different from that of 
tungsten or molybdenum. While tantalum can be heated to as high, or higher, a tem- 
perature than molybdenum, the adsorbed gases are not as readily desorbed. 

_Q 

Hydrogen: At a hydrogen pressure of 5x10 torr, the desorption curves of Hg-Ta 
were similar to those obtained on Mo or W, except that the temperature level was con- 
siderably higher. A typical trace is shown in sketch (m). The temperature of the maxi- 
imum desorption rate is plotted in figure 5(b) for several such traces. The general 
desorption energy level from this data is 2. 9 to 3. 1 electron volts assuming second- 
order desorption. 

A plot of the buildup of surface concentration with adsorption time at an H 0 pressure 
of ~lxl0" torr is shown in sketch (n). A maximum coverage of 1. 56xl0 10 H-atoms 
per square centimeter under these conditions was measured. This number is in the 
range that might be considered monolayer coverage, though it is higher than anticipated. 
However, when hydrogen pressures up to 10” 5 torr were introduced, the subsequent 
desorption curves yielded amounts of desorbed H 2 much in excess of this "monolayer 
range" value. In fact, desorbed H 2 equivalent to nearly 50 times this value was mea- 
sured. Even this was not necessarily a limiting value, but further pressurization was not 
done. It seems that this excessive amount of hydrogen must be penetrating into the in- 
terior of the tantalum metal. However, the temperature of maximum desorption rate 
showed a decreasing trend with "coverage" as shown in sketch (o). This trend seems 
peculiar in that, if hydrogen were penetrating into the metal lattice, the subsequent time 
required for the H 2 to reach the surface before desorbing would be expected to show up as 
as a higher maximum temperature level in the desorption trace. A desorption energy of 
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2. 82 electron volts is reported for this system by Krakowski (ref. 16). He also suggests 
Hg solubility in bulk tantalum as a possible explanation of some of his observed behavior. 

Other gases: It was of interest to see if the adsorption characteristics of deuterium 
were the same as Hg on tantalum. The temperature range of desorption and the temper- 
ature of maximum desorption rate for a number of desorption traces indicated that Dg 

behavior was the same as H 0 with tantalum. T values of 1250±100 K were obtained 

for over a dozen desorption traces, indicating the Dg-Ta desorption energy level to be 
around 2. 8 to 3. 2 electron volts. 

Carbon monoxide desorption from tantalum started at a temperature level of around 
1900 K, but never reached a maximum desorption rate for wire temperatures up to 
about 2800 K. 

No N 9 desorption was observed from Ta, even though adsorption pressure levels as 
-4 

high as 10 torr were used. 

The fact that CO desorption was never complete, and no nitrogen desorption observed 
indicates that the tantalum surface was probably never very clean. Then, the hydrogen 
and deuterium desorption cited previously is probably from some Ta-X surface. 


Site-Energy Determination 

The experimentally determined curves of surface coverage ratio ®/ Q m against 
temperature are shown in sketches (a) and (p). 

These experimental curves for completely covered monolayers for the systems 
CO-W, CO-Mo, Hg-W, Ng-Mo, and Ng-W were used to calculate energy-area relations 
as previously outlined. The results of these calculations are shown in figure 7. 

Also shown in these figures are the F(E) curves which were determined directly 
from the graphical plots (sketches (a) and (p)). It is apparent that the graphical approach 
outlined yields a reasonable approximation to the computer -determined distribution, 
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except for the case of the ^-W data. For these cases, at least, the graphical approach 
yields a slightly narrower energy spread than the computer solution. As was mentioned 
previously, the effect of plateaus in the F(E) against E curves may not be readily 
apparent in the desorption traces. It is not certain that the inflections in the energy 
distribution curve for the Ng-W system are representative of the gas -surface interaction 
energy entirely for that system, or are affected by the sensitivity of the particular 
numerical calculation. The resulting isotherm calculations will not be greatly in- 
fluenced by small variations in the site -energy distribution vurves. 

This experiment alone is not capable of distinguishing between variations in E(6) 
due to interactions between adsorbed particles and the variations due to surface hetero- 
geneity. In order to see if there were any similarities in desorption behavior that might 
be attributed to surface variation alone, the F(E) curves of figure 7 were plotted with the 
energy scale normalized to the E -value at half coverage (sketch (q)). The shapes of 
these curves are different for molybdenum and tungsten, but do have similarity between 
different gases on the same surface, which is at least suggestive of a given surface - 
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patch distribution. The Ng-W data again are not consistent in this comparison. 

The general spread in desorption energies for these cases is seen to be from about 
30 percent below to about 40 percent above the energy at half coverage. 


Isotherm Determination 

Sticking coefficient . - For conditions of negligible desorption rate, the. sticking co- 
efficient is defined as follows: 


Adsorption rate = juS (2) 

The time derivative of the surface coverage (for instance, the slope of the curve in 
sketch ( l ) or (o)) can give values of S. If, at the values of coverage and temperature 
involved, the desorption rate is not negligible, then the determination of S from ex- 
perimental data must include the desorption rate term in the relation; that is, 

Adsorption rate = fiS - v (42) 

Since v is a relatively sensitive function of T_ and E, as well as coverage, it is 
obviously less accurate to determine a value of sticking coefficient under these condi- 
tions. 

Generally, the sticking coefficient is calculated through equation (2) even though 
the conditions of negligible desorption rate may not be met. This is done because the 
experiment frequently simply does not yield the attendant values of v to use. This 
leads to "pseudosticking" coefficients which drop to zero at coverages of less than 
unity. 

Sticking coefficients determined from data presented herein are shown in sketch (r). 
The solid lines are plots of the relation S = S Q (1 - 9), with S Q = 0. 37 for ^-W, 

0. 17 for Hg-Ta, and 0. 22 for Hg-W. (Comparative values of S from the literature are 
also given in table n). 

There is considerable variation in values of S in the literature for a given system 
(see ref. 17, for example). Some probable reasons for the variations are also cited 
therein. A wide variation in sticking coefficient with crystal face has been reported 
for the Hg-W system (ref. 18). Presence of impurity gases can cause error in the 
measurements. In view of the relative uncertainty of sticking coefficient values, the 
isotherm calculations presented here were made using only the clean surface value of 
sticking coefficient S Q . Variation with coverage is assumed to follow equation (6). The 
trend of values of S shown in sketch (r) as well as other reported values (ref. 17) 
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indicates this to be a reasonable assumption. The assumption is the same as that used 
in calculating the generalized isotherm plots of figure 1. 

Calculated isothe rms. - Isotherm data calculated for the same systems are shown 

in figure 8. These isotherms were determined using the computer-determined f(E) re 

lations and equation (24). Other data necessary to calculate these plots are included in 

12 

table I. Pressure range selected for these plots was from 10 to 1 torr with the 
appropriate temperature range to yield fractional coverages from 0.02 to 0.98. 

Lines of constant total coverage 0 did plot as straight lines on the log p against 
1/T plots. Thus, a constant value of desorption energy E(0) can be used to express 
the total desorption rate from the heterogeneous surface at any given total coverage 
value. The values of E(0) determined from the slopes of these lines are shown in 
figure 7 also for comparison. Curves for E(0) against 0 are plotted in figure 9. 

The graphical method of obtaining the F(E) against E relation produces values 
that may be used to describe the distribution relation and to calculate adsorption iso- 
therms. It appears that, if the true F(E) is a reasonably continuous function of energy 
the resulting calculations using the graphical F(E) relation will be a good approxima- 
tion to the presumably more correct computer values. 


CONCLUDING REMARKS 

Equations and curves describing the adsorption of gas on metal substrates under 
equilibrium conditions, and the desorption of gas from the surface under a steady tem- 
perature increase of the surface temperature have been presented. 

Results of an experimental investigation using the flash -filament desorption tech- 
nique are presented. Gas -metal surface combinations for which quantitative experi- 
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mental data are presented are: hydrogen-tungsten, carbon monoxide -tungsten, 
hydrogen-molybdenum, carbon monoxide-molybdenum, nitrogen-molybdenum, and 
hydrogen-tantalum. Some observations on the behavior of other gas -metal combinations 
were also made. 

The experimental data were treated to yield energy-distribution functions for the 

gas -metal combinations. A simplified graphical method of approximating the energy- 

distribution function was presented. This method gave values which can be used to 

obtain adsorption isotherms in good agreement with the computer -determined values. 

Adsorption isotherms were calculated using the experimentally determined energy- 

distribution relations. These isotherm plots along with the monolayer coverage value 

or express the expected equilibrium amount of gas adsorbed over a range of pressures 
m .10 

from 10 to 1 torr, surface coverages from 0.02 to 0. 98 monolayers, and for the 
appropriate temperature range in each case. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 6, 1968, 

129-01-02-07-22. 


27 


I 



APPENDIX 


SYMBOLS 


a 


r d 

F(E),f(E) 

f(jtt,T s ) 

h 

k 

L 

m 

N 

n 



^exh 


S 

S c 

T 

t 

V 




surface area, cm^ 
transmission area, cm^ 
wire diameter, cm 

desorption energy, eV (1 eV = 1. 6x10"*^ J) 

net flow out of test chamber, particles/sec 

combined flow into test chamber, particles/sec 

desorption rate, particles/sec 

adsorption rate, particles/sec 

site -energy distribution functions 

function of fj. and T_ 

s 

Planck’s constant, J-sec 
Boltzmann constant, eV/K 
wire length, cm 
molecular weight 
number of molecules 
order of reaction 
pressure 

chamber pressure, torr 

equilibrium pressure in test chamber 

pressure of exhaust chamber 

sticking coefficient 

sticking coefficient at zero coverage 

temperature 

time 

volume 

rate of temperature increase of surface, dT/dt, K/sec 
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© overall fractional monolayer coverage of heterogeneous surface 
9 fractional monolayer coverage, CT/cr m , of homogeneous patch 

9 

fx arrival rate, particles/(cm )(sec) 

v desorption rate, particles/(cm )(sec) 

2 

a surface coverage, particles/cm 

2 

a m monolayer surface coverage, particles/cm 

2 

0q initial surface coverage, particles/cm 

t time constant in eq. (20), sec 

T on constant in n m -order desorption equation, sec 

t i constant in first-order desorption equation, sec 

t o 2 constant in second-order desorption equation, sec 

Subscripts: 

g gas 

i i** 1 component of area 

max maximum 

s surface 

1,2 initial and final states of process 
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TABLE I. - SUMMARY OF PARAMETERS 


System 

Maximum coverage, 

°max ? 

-2 

cm 

Clean surface 
sticking coef- 
ficient, 

S o 

Temperature of 
maximum desorp- 
tion rate, 

T 

max’ 

K 

Desorption energy 
associated with 
maximum rate, 
E, 
eV 

Reaction order, 
n 

h 2 -w 

1. OxlO 15 H -atoms 

0.22 

740 to 860 

1.7 to 2.0 

2 

n 2 -w 

6. 44X10 14 N-atoms 

.037 

1500 to 1700 

3. 7 to 4. 1 

2 

co-w 

5. 3xl0 44 CO-molecules 

a . 5 

1370 to 1430 

3. 3 to 3. 4 

1 

N 2 -Mo 

1. 16X10 15 N-atoms 

b .7 

1180 to 1220 

3.0 

2 

CO-Mo 

4. 83xl0* 4 CO-molecules 

c 1.0 

1220 to 1270 

3. 1 to 3. 3 

1 

H 2 -Ta 

**1. 58xl0 4 ® H-atoms 

.17 

1170 to 1280 

2.9 to 3.1 

2 

H 2 -Mo 

e 2. 65xl0 14 H-atoms 

.24 

570 to 700 

1.35 to 1.55 

2 


‘‘Ref. 19. 
b Ref. 20. 

c Estimated in lieu of experimental value. 
^Not maximum coverage, see text. 
e Not believed maximum attainable. 
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TABLE II. - DESORPTION ENERGY, STICKING COEFFICIENT, AND 
MONOLAYER COVERAGE FOR SOME GAS-METAL SYSTEMS 
[All values as cited in ref. 2, but only those of interest above 300° are cited. ] 


System 

Desorption energy, 
E, 
eV 

Clean surface sticking 
coefficient, 

S o 

Monolayer surface 
coverage, 
cr , 

nr 2 

particles/cm 

H 2 -w 

1.35, 2.9, 1. 48±0. 07, 
1.10, 1.26, 0.87 to 2.0 

0.1, 0.2 
0. 1 (100 plane) 

0. 3 (211 plane) 
-0.01 (111 plane) 
<10" 3 (110 plane) 

7.5X10 14 

n 2 -w 

3.5 

0. 3 to 0. 4 

3xl0 14 to 10xl0 14 

co-w 

1 

a v 0.87 to 1. 26 a 
2.26 to 2.74 a 
/3 2 , 3.02 to 3. 35 a 
j3g, 3.28 to 4. 35 a 

0.2, 0.49, 0.50, 0.48 

4. 5xl0 14 to 9.5X10 14 

H 2 Mo 

i 

j 

0.35 

8.2X10 14 

N 2 -Mo 

1 

i 

0.6, 0.7 

. i 

4.4X10 14 


designations a and (3 are as used in ref. 2. 



Desorption energy, E, eV Coverage fraction, 0 Coverage fraction, 





Figure 1. - Variation of equilibrium surface coverage with arrival rate and surface temperature 
for surfaces of constant desorption energy. First- or second-order desorption. 







Figure 3. - 
tion of coi 


(a -2) Temperature rise rate, 100 K per second. 



(a-3) Temperature rise rate, 1000 K per second. 

(a) First-order desorption, n * 1; t o1 = 10" 13 

ation of surface coverage with surface temperature during desorp 
it energy surface at various temperature rise rates. 










O Hydrogen-molybdenum 
A Hydrogen -tungsten 
D Hydrogen -tungsten (simulta- 
neous carbon monoxide) 

O Nitrogen -tungsten 
^ Nitrogen-molybdenum 
V Hydrogen -tantalum 

Range of oxygen -tungsten data 



Temperature of maximum desorption rate. T max , K 


(a) First-order desorption, n - 1. 



200 600 1000 1400 1800 2200 2600 3000 

Temperature of maximum desorption rate, T max , K 


(b) Second-order desorption, n = 2. 

Figure 5. -Temperature at which maximum desorption rate occurs for constant 
energy surface. 
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(b) Tungsten. 

Figure 6. - Desorption curves for maximum coverage 
cases. 
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Site-energy function, F(E) 




(a) Carbon monoxide -tungsten. 


(b) Carbon monoxide-molybdenum. 




(c) Hydrogen -tungsten. 


(d) Nitrogen -molybdenum. 
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Desorption energy, E, eV 


(e) Nitrogen -tungsten. 

Figure 7. - Energy distribution function variation with energy. 







Temperature, T, K 
(e) Carbon monoxide-molybdenum. 
Figure 8. - Concluded. 



Figure 9. - Desorption energies determined from slopes of constant 
coverage lines of figure 8. 
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